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Method and arrangement for controlling injection of reducing agent 
FIELD OF THE INVENTION, AND STATE OF THE ART 

10 The present invention relates to a method and a device for 
controlling the injection of reducing agent upstream from a 
catalyst in an exhaust line from a combustion engine. The 
invention also relates to a computer program comprising program 
codes for implementation of said method, a computer program 

15 product comprising a medium which is readable by an electronic 
control unit and has stored on it a computer program intended to 
cause an electronic unit to implement said method, and an 
electronic control unit. 

20 To meet prevailing exhaust cleaning requirements, today's motor 
vehicles are usually provided with a catalyst in the exhaust line 
for effecting catalytic conversion of environmentally harmful 
components of exhaust gases to environmentally less harmful 
substances. One method employed for achieving effective 

25 catalytic conversion is based on injection of a reducing agent in 
the exhaust gases upstream from the catalyst. A reducing 
substance forming part of or formed by the reducing agent is 
carried by the exhaust gases into the catalyst, in which it is 
adsorbed on active seats in the catalyst, thereby causing storage 

30 (accumulation) of the reducing substance in the catalyst. The 
stored reducing substance may either desorb, i.e. become 
detached from the active seats, or react with an exhaust gas 
substance so as to convert said exhaust gas substance to a 
harmless substance. Such a reduction catalyst may for example 

35 be of SCR type (SCR = Selective Catalytic Reduction). This type 
of catalyst is hereinafter called SCR catalyst. An SCR catalyst 
selectively reduces NO x in exhaust gases but not the oxygen in 
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exhaust gases. In the case of an SCR catalyst, a reducing agent 
in the form of urea or ammonia is usually injected in the exhaust 
gases upstream from the catalyst. Injection of urea into exhaust 
gases results in the formation of ammonia and it is this ammonia 
5 that constitutes the reducing substance which assists the 
catalytic conversion in the SCR catalyst. The ammonia is 
accumulated in the catalyst by being adsorbed on active seats in 
the catalyst, and NO x present in the exhaust gases is converted 
to nitrogen and water by being brought into contact with ammonia 
10 accumulated on active seats in the catalyst. 

When a reduction catalyst is used in combination with 
proportioning of reducing agent, it is important to control the 
injection of the reducing agent so that desired conversion of the 

15 exhaust gas substance concerned is achieved without allowing 
excessive amounts of unconsumed reducing substance to 
accompany the exhaust gases leaving the catalyst and thereby 
reach the environment. The consumption of reducing agent and 
the relating costs can thereby be minimised while at the same 

20 time preventing or at least minimising undesired release of the 
reducing substance to the environment. The fact, for example, 
that ammonia is a reducing substance which is toxic and 
malodorous makes it desirable to minimise to the utmost possible 
extent the release of this ammonia to the environment while 

25 maintaining the injection of sufficient reducing agent, e.g. in the 
form of urea, for achieving the required conversion of the 
exhaust gas substance concerned, in this case NO x . 

When a catalyst of the type here concerned is used, the injection 
30 of reducing agent is at the present time usually controlled on the 
basis of two-dimensional tables, each of which is related to a 
certain exhaust gas temperature. Each table comprises values 
for the amount of reducing agent to be injected in the exhaust 
line at the respective load and speed of the combustion engine 
35 concerned. Each table thus has engine load on a first axis and 
engine speed on a second axis. These tables are compiled 
empirically by charting the behaviour of the respective 
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combustion engine and catalyst on the basis of performance in 
steady-state conditions. Such charting has to be done for each 
specific combination of combustion engine and catalyst, and this 
previously known solution also involves extensive trial runs and 
5 settings for each specific combination of engine type and catalyst 
type. Any changes in, for example, engine emissions entail 
having to make extensive and detailed changes to all the values 
in the tables. As the tables give proportioning adapted to steady- 
state conditions, dynamic compensations have to be applied in 
10 cases where operating conditions change quickly. A further 
disadvantage of the known solution is that, despite dynamic 
compensations, operating conditions for which the tables do not 
cater can occur, which may lead to incorrect proportioning of 
reducing agent. 

15 

OBJECT OF THE INVENTION 

The object of the present invention is to provide a method which 
enables easy and reliable control of the injection of reducing 
20 agent upstream from a catalyst in an exhaust line from a 

combustion engine so that amounts of reducing agent which are 
appropriate to prevailing operating conditions can be injected. 

SUMMARY OF THE INVENTION 

25 

According to the invention, said object is achieved by means of a 
method exhibiting the features indicated in claim 1. 

The solution according to the invention involves: 
30 - an accumulation actual value representative of current 

accumulation in the catalyst of a reducing substance forming part 
of or formed by the reducing agent being calculated on the basis 
of information from a computation model which, taking into 
account the reactions in the catalyst under prevailing operating 
35 conditions, continuously determines the current state of the 
catalyst, e.g. the accumulation of the reducing substance in 
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different parts of the catalyst and the conversion of exhaust gas 
substance taking place in different parts of the catalyst, 

- an accumulation setpoint value being calculated on the basis of 
an emission setpoint value and information from said 

5 computation model, whereby the emission setpoint value is 

representative of a desired content, in the exhaust gases leaving 
the catalyst, of an exhaust gas substance which, as the exhaust 
gases pass through the catalyst, is at least partly removed from 
the exhaust gases by the action of the reducing substance or is 

10 formed by the action of the reducing substance, and the 

accumulation setpoint value is representative of the accumulation 
of the reducing substance which is required in the catalyst under 
prevailing operating conditions for substantially achieving the 
emission setpoint value, 

15 - the accumulation actual value being compared with the 
accumulation setpoint value, and 

- the injection of reducing agent in the exhaust line being 
controlled on the basis of the result of the comparison between 
the accumulation actual value and the accumulation setpoint 

20 value. 

The solution according to the invention makes it possible easily 
and reliably to control the injection of reducing agent upstream 
from a catalyst in an exhaust line from a combustion engine so 

25 that suitable amounts of reducing agent can be injected, based 
on prevailing operating conditions, for achieving a desired 
conversion of the exhaust gas substance concerned. Taking the 
current reducing substance accumulation in the catalyst into 
account makes it possible to control the proportioning of reducing 

30 agent with good accuracy so that the desired conversion of the 
exhaust gas substance is accomplished while at the same time 
the amount of unconsumed reducing substance which 
accompanies the exhaust gases leaving the catalyst and thereby 
reaches the environment is limited. It thus becomes possible to 

35 minimise the consumption of reducing agent and the relating 
costs while at the same time preventing or at least minimising 
undesired release of reducing agent to the environment. The 
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computation model used need only be adapted to the catalyst 
and is therefore independent of the design of the combustion 
engine connected to the catalyst. Using a computation model 
which is only tied to the catalyst and not to the combustion 
5 engine results in a very flexible method which, when adapted to a 
specific catalyst, can be used for that catalyst together with any 
desired combustion engine. It is also possible to use a scalable 
computation model which can easily be adapted to the 
dimensions of the catalyst. 

10 

According to an embodiment of the method according to the 
invention, a limitation factor is calculated, which limitation factor 
has a value based on an estimate of the current risk that the 
reducing substance content of the exhaust gases leaving the 

15 catalyst might exceed a predetermined limit value, whereby this 
limitation factor is taken into account in calculating the 
accumulation setpoint value in such a way that the accumulation 
setpoint value decreases in response to increasing risk that the 
reducing substance content of the exhaust gases leaving the 

20 catalyst might exceed the predetermined limit value. It thus 
becomes possible to eliminate or at least minimise the risk that 
undesired amounts of unconsumed reducing substance might 
accompany the exhaust gases leaving the catalyst and thereby 
reach the environment. 

25 

A further embodiment of the method according to the invention is 
characterised by: 

- an emission actual value being established by calculation or 
measurement, which emission actual value is representative of 

30 the current content of the exhaust gas substance in the exhaust 
gases leaving the catalyst, 

- the emission actual value being compared with the emission 
setpoint value, and 

- the accumulation setpoint value being calculated on the basis 
35 of information from said computation model and the conformity 

between the emission actual value and the emission setpoint 
value. 
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The fact of thus causing the conformity between the emission 
actual value and the emission setpoint value to affect the 
accumulation setpoint value makes it possible to calculate the 
accumulation setpoint value by an approximation method, since 
5 the influence on the accumulation setpoint value of the 

conformity between the emission actual value and the emission 
setpoint value makes it possible to compensate for deviations 
arising between desired and actual emissions of the exhaust gas 
substance. The calculated accumulation setpoint value thus 
10 need not correspond exactly to the reducing substance 

accumulation required for achieving the emission setpoint value. 

Further embodiments of the method according to the invention 
are indicated by the dependent claims and the description set out 
15 below. 

The invention also relates to a device according to claim 21 for 
utilising the method according to the invention. 

20 The invention also relates to a data program according to claim 
23 which is loadable directly to the internal memory of a 
computer and which comprises program codes for 
implementation of the method according to the invention. 

25 The invention also relates to a computer program product 

according to claim 24 comprising a medium which is readable by 
an electronic control unit and which has stored on it a computer 
program intended to cause an electronic control unit to 
implement the method according to the invention. 

30 

The invention also relates to an electronic control unit according 
to claim 25. 

BRIEF DESCRIPTION OF THE DRAWINGS 



35 
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The invention is described in more detail below on the basis of 
examples of embodiments with reference to the attached 
drawings, which are as follows: 

5 Fig. 1 a schematic diagram of a combustion engine with a 
relating catalyst and a device according to the 
invention, 



Fig. 2 a block diagram illustrating a first embodiment of the 
10 device according to the invention, 

Fig. 3 a block diagram illustrating a second embodiment of 
the device according to the invention, 

15 Fig. 4 a block diagram illustrating an electronic control unit 
for implementation of the method according to the 
invention, 

Fig. 5 a flowchart illustrating a method according to a first 
20 embodiment of the present invention, 

Fig. 6 a flowchart illustrating a method according to a second 
embodiment of the present invention, 

25 Fig. 7 a flowchart illustrating a method according to a third 
embodiment of the present invention, and 



30 



Fig. 8 a flowchart illustrating a method according to a fourth 
embodiment of the present invention, and 

Fig. 9 a flowchart illustrating a method according to a fifth 
embodiment of the present invention. 



DETAILED DESCRIPTION OF EMBODIMENTS 

35 

Fig. 1 depicts schematically a combustion engine 1 provided with 
a device 10 according to the invention. The exhaust gases 
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leaving the combustion engine 1 pass through an exhaust line 2 
and reach the environment via an exhaust outlet 3. A reduction 
catalyst 4, preferably in the form of an SCR catalyst, is arranged 
in the exhaust line 2. The exhaust gases from the combustion 

5 engine 1 are thus caused to pass through this catalyst 4 before 
reaching the environment via the exhaust outlet 3. An injection 
point 5 for reducing agent is situated in the exhaust line 2 
upstream from the catalyst 4. The injection of reducing agent 
takes place via an injection device comprising one or more 

0 injection means 6 arranged in the exhaust line, in the form of 
injection nozzles or the like, and a reducing agent storage tank 7 
connected thereto. The injection device further comprises a 
proportioning unit 8, which itself comprises proportioning means 
and a regulating means, e.g. a control valve or the like, arranged 

5 to regulate the supply of reducing agent to said injection means 
6. The proportioning unit 8 is controlled by a device 10 
according to the present invention which determines the amount 
of reducing agent to be injected in the exhaust gases. 

!0 Fig. 2 illustrates a first embodiment of a device 10 according to 
the invention. The device 10 comprises a first calculation means 
20 adapted to determining continuously, by use of a computation 
model, the current state of the catalyst, taking into account the 
expected reactions in the catalyst under prevailing operating 

>5 conditions. The computation model takes into account the 

expected exothermal and endothermal reactions in the catalyst 4 
under prevailing operating conditions and calculates inter alia the 
accumulation of the reducing substance concerned in different 
parts of the catalyst, and the conversion of the exhaust gas 

30 substance concerned which takes place in different parts of the 
catalyst. The computation model may be of any desired design 
provided that it produces with the desired accuracy a correct 
value for the accumulation of the reducing substance and the 
conversion of the exhaust gas substance in the catalyst. An 

35 example of a computation model suitable in the context is 
described below. 
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The device 10 further comprises a second calculation means 30 
adapted to calculating, on the basis of information from the first 
computation means 20, an accumulation actual value A 1 
representative of the current accumulation in the catalyst 4 of the 
5 reducing substance which forms part of or is formed by the 
reducing agent. The second computation means 30 is also 
adapted to calculating an accumulation setpoint value A2 on the 
basis of an emission setpoint value E2 and information from said 
computation model. This emission setpoint value E2 is 

10 representative of a desired content, in the exhaust gases leaving 
the catalyst 4, of an exhaust gas substance which, as the 
exhaust gases pass through the catalyst, is at least partly 
removed from the exhaust gases by the action of the reducing 
substance or is formed by the action of the reducing substance. 

15 The emission setpoint value E2 provides a measure of the 

maximum permissible emission of the exhaust gas substance and 
may inter alia be controlled by statutory exhaust emission 
requirements. The emission setpoint value E2 has in the 
simplest case a constant value but is preferably determined by 

20 calculation in a computation means 38 in a conventional manner 
on the basis of prevailing operating conditions, e.g. the speed 
and load of the combustion engine. The accumulation setpoint 
value A2 is representative of the reducing substance 
accumulation required in the catalyst under prevailing operating 

25 conditions for achieving or at least substantially achieving the 
emission setpoint value E2. The first computation means 20 and 
the second computation means 30 are with advantage integrated 
in a common computer unit but may, if considered appropriate, 
take the form of separate units connected to one another. 

30 

The device 10 further comprises a comparator 32 adapted to 
receiving the accumulation actual value A1 and the accumulation 
setpoint value A2 from the second computation means 30. The 
comparator 32 is adapted to comparing the accumulation actual 
35 value A1 and the accumulation setpoint value A2 and to emitting 
a signal S1 which depends on the conformity, e.g. the difference, 
between this actual value A1 and setpoint value A2. The device 
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10 also comprises regulating means 34, 36 for controlling the 
injection of reducing agent on the basis of said signal S1 from 
the comparator 32. Said regulating means comprise with 
advantage a regulator 34, preferably in the form of a PI regulator 
5 or PID regulator, adapted to receiving the signal S1 from the 
comparator 32. On the basis of that signal S1, the regulator 34 
emits a control signal S2 to a control means 36 connected to the 
proportioning unit 8 of the injection device and adapted to 
controlling this proportioning unit on the basis of said control 
10 signal S2 so that an amount of reducing agent appropriate to 
prevailing conditions is injected in the exhaust line 2. 

The control signal S2 from the regulator 34 may also be led to 
the first computation means 20 in order to be used in the 

15 computation model, as a measure of the amount of reducing 
agent injected, when generating information for the next 
calculation of accumulation actual value A1 and accumulation 
setpoint value A2. Alternatively, the first computation means 20 
may be adapted to receiving from the control means 36 a value 

20 representing the amount of reducing agent injected, as illustrated 
in Fig. 3. 

Fig. 3 illustrates a second embodiment of a device 10 according 
to the invention. The components which appear both in the 

25 embodiment described above with reference to Fig. 2 and in the 
embodiment according to Fig. 3 are denoted by the same 
reference notations. In the embodiment according to Fig. 3, the 
device 10 likewise comprises means for determining by 
calculation or measurement an emission actual value E1 

30 representative of the current exhaust gas substance content of 
the exhaust gases leaving the catalyst 4. In cases where the 
emission actual value E1 is determined by calculation, this is 
preferably done in the first computation means 20 by using the 
abovementioned computation model or in a separate computation 

35 means on the basis of information from that computation model. 
In cases where the emission actual value E1 is determined by 
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measurement, this is done by a measuring means arranged in the 
exhaust line downstream from the catalyst 4. 

In the embodiment according to Fig. 3, the second computation 
5 means 30 is adapted to calculating the accumulation setpoint 
value A2 on the basis of information from the first computation 
means 20 and the conformity between the emission actual value 
E1 and the emission setpoint value E2. In this case the device 
10 comprises a comparator 42 adapted to receiving the emission 

10 actual value E1 and the emission setpoint value E2. The 

comparator 42 is adapted to comparing the emission actual value 
E1 and the emission setpoint value E2 and to emitting a signal 
S3 which depends on the conformity, e.g. the difference, 
between this actual value E1 and setpoint value E2. The device 

15 10 comprises with advantage a regulator 44, preferably in the 
form of a PI regulator or PID regulator, which is adapted to 
receiving the signal S3 from the comparator 42. On the basis of 
that signal S3 the regulator 44 emits a control signal f S p to the 
second computation means 30, which then determines the 

20 accumulation setpoint value A2 on the basis of that control signal 
f S p and the information from the first computation means 20 so 
that the accumulation setpoint value A2 is to the necessary 
extent compensated for undesired deviations between the 
emission actual value E1 and the emission setpoint value E2. 

25 

In the embodiment according to Fig. 3, regulation is applied in 
the form of so-called cascade regulation whereby the second 
computation means 30, the comparator 32 and the regulator 34 
constitute an inner regulating circuit 12a which sets a necessary 

30 reducing substance accumulation in the catalyst by acting upon 
the proportioning of the reducing agent, and whereby the 
comparator 42 and the regulator 44 constitute an outer regulating 
circuit 12b which sets a desired emission of the exhaust gas 
substance by acting, via the inner regulating circuit 12a, upon the 

35 level of the necessary accumulation of the reducing substance in 
the catalyst. 
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Program codes for implementation of the method according to the 
invention are preferably adapted to forming part of a computer 
program directly loadable to the internal memory of a computer, 
e.g. to the internal memory of the aforesaid computer unit. Such 
5 a computer program is provided with advantage via a computer 
program product comprising a storage medium which is readable 
by a computer and on which the computer program is stored. 
Said storage medium takes the form, for example, of an optical 
storage medium in the form of a CD-ROM disc, a DVD disc etc., 

10 or of a magnetic storage medium in the form of a hard disc, a 
diskette, a cassette tape etc. Fig. 4 illustrates an electronic 
control unit 50 comprising a means 51, preferably a central 
processor unit (CPU), for execution of software, which 
communicates via a databus 52 with a memory 53, e.g. of the 

15 RAM type. The control unit 50 also comprises at least one 
storage medium 54, e.g. in the form of a memory of the ROM, 
PROM, EPROM or EEPROM type or a Flash memory, which the 
execution means 51 communicates with via the databus 52. A 
computer program comprising program codes for implementation 

20 of the method according to the invention is stored in the storage 
medium 54. 

The reducing agent preferably takes the form of urea (CO(NH 2 )2) 
but may also take the form of, for example, ammonia (NH 3 ) or 

25 hydrocarbon (fuel). In the ensuing description it is assumed that 
the catalyst 4 is an SCR catalyst and that a reducing agent in the 
form of urea or ammonia is used. It should be noted, however, 
that the solution according to the invention is not limited to this 
type of reduction catalyst and this type of reducing agent. 

30 Injecting urea into exhaust gases results in the formation of 
ammonia and it is this ammonia that constitutes the reducing 
substance which assists the catalytic conversion in the SCR 
catalyst. In this case the exhaust gas substance concerned 
takes the form of NO x . 

35 

The aforesaid computation model is with advantage designed to 
use the following parameters as input values: 
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a) The exhaust gas temperature P1 upstream from the catalyst 4. 
This temperature may be determined by temperature sensor or 
in any of the conventional ways of calculating it. 

5 

b) The NOx concentration P2 in the exhaust gases upstream from 
the catalyst 4. This concentration may be determined by 
sensor but is determined with advantage in any of the 
conventional ways of calculating it, e.g. on the basis of the 

10 combustion engine's load, speed, injection angle, i.e. the 

angle of the combustion engine's crankshaft at the time of fuel 
injection into the engine cylinders, and, where applicable, the 
EGR content (EGR = Exhaust Gas Recirculation), i.e. the 
content of exhaust gases led back to the engine. 

15 

c) The exhaust gas mass flow P3 through the catalyst 4. This 
exhaust gas mass flow may be determined by mass flow 
sensor but is determined with advantage in any of the 
conventional ways for calculating it, e.g. on the basis of the 

20 combustion engine's load and speed. 

d) The amount P4 of reducing agent injected in the exhaust 
gases. The value for the amount of reducing agent injected is 
obtained with advantage from the control means 36 or from 

25 the regulator 34. 

The computation model may also use as input value the 0 2 
concentration in the exhaust gases upstream from the SCR 
catalyst and/or the ambient temperature. The 0 2 concentration 
30 may be determined by, for example, lambda sensor but is 

determined with advantage in any of the conventional ways for 
calculating it, e.g. on the basis of the combustion engine's load, 
speed and, where applicable, EGR content. 

35 A design of computation model is described below which it is 

advantageous to use for obtaining information for calculating the 
aforesaid accumulation actual value A1, accumulation setpoint 
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value A2 and, where applicable, emission actual value E1 which 
pertain to a method and a device according to the present 
invention. 

5 In an SCR catalyst, nitrogen oxide NO x reacts with ammonia and 
is reduced to nitrogen gas. NO x is the harmful exhaust gas 
substance intended to be removed from the exhaust gases, and 
ammonia is the reducing substance used for achieving this. 
Ammonia or urea (which converts to ammonia) is injected in the 

10 exhaust gases upstream from the SCR catalyst. The 

computation model is used to determine how much NO x is 
converted in the SCR catalyst and how much unconsumed 
ammonia leaves the SCR catalyst. The computation model also 
calculates continuously how the temperature varies across the 

15 catalyst and how much ammonia is accumulated in different parts 
of the catalyst. This involves the computation model being 
continuously supplied with information about the magnitude of 
the gas flow through the catalyst and the temperature and 
composition of the gas flowing into the catalyst. 

20 

A number of reactions take place in the SCR catalyst. Ammonia 
is adsorbed on active seats in the catalyst, resulting in 
accumulation of ammonia in the catalyst. The accumulated 
ammonia may either desorb, i.e. become detached from the 
25 active seats, or react with NO x . At high temperatures, oxidation 
of ammonia with oxygen also takes place to a certain extent. 
The factor which decides how much NO x is converted in the 
catalyst is the reaction velocities r, of the various reactions. The 
reactions and their relating reaction velocities are as follows: 

30 

1) S + NH 3 S.NH 3 r, = k^naOv 

2) S.NH 3 -^ S + NH 3 r 2 = k 2 9 NH3 

35 3) 4S.NH 3 + 4NO + Q 2 -*- 4S + 4N 2 + 6H 2 0 r 3 = k 3 c N oO NH3 



4) 4S.NH 3 + 5Q 2 4S + 6H 2 0 + 4NO r 4 = k 4 Co 2 e NH3 
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where k, is the velocity constant for reaction /, c, is the 
concentration of substance /', 9 V is the proportion of vacant seats, 
and 9nh3 is the proportion of seats occupied by ammonia. The 
5 reaction velocities r, are temperature-dependent in accordance 
with the Arrhenius equation: 

ki = k 0J e «t 

10 where k 0J is constant for reaction /, E AJ is the activation energy 
for reaction /, R is the general gas constant, and T is the 
temperature. 

To determine the composition of the exhaust gases after the 

15 catalyst and how much ammonia is accumulated in different parts 
of the SCR catalyst, a number of material balances are solved 
according to the computation model. As the SCR catalyst has a 
monolithic structure, the gas flows through small channels 
whereby the walls between the channels contain the active 

20 catalyst material. The catalyst is modelled by regarding the flow 
through a channel as divided into a number of segments. The 
material balances are solved successively from the segment at 
the catalyst's inlet to the segment at the catalyst's outlet. From 
the flow through the channel, MO x and ammonia are transported 

25 up to the wall of the duct, where these substances react. To take 
into account the effect of the velocity at which the substances 
are transported up to the channel wall and into the channel wall, 
the channel wall is likewise divided into a number of segments. 
As all the material balances in the wall segments within each 

30 channel segment are inter-related, they have to be solved 

together in an equation system. According to the computation 
model, the following material balances are set up: 

f uh O'a-i - y ijc ) - r a>0 (c itkfi - c w ) = o 

35 

ru >B _, (<?,,*,„_, -c a ,„) -r,. M (c fji>n -c,. t n+1 ) + Yj v u r j*.n w k. n =o for /7 > 1 
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where F tot is the total molar flow, y lik and c, ik are the molar 
proportion and the concentration respectively of substance / in 
channel segment k, r ao and T iJcn are the respective coefficients 

5 for transport of substance / from the gas flow to the first wall 
segment and between wall segments n and n+1 in channel 
segment k, v u are stoichiometric coefficients for substance / in 

reaction j, r jikn is the reaction velocity for reaction j in channel 
segment k and wall segment n, and w kw „ is the mass of active 
10 catalyst material in channel segment k and wall segment n. The 
ammonia accumulation in channel segment k and wall segment n 
is then arrived at by the material balance: 

N de ™*** = y v r 

15 

where N c is the number of active seats per unit mass of catalyst. 

To determine the temperature across the SCR catalyst, a heat 
balance for the gas and a heat balance for the catalyst are 
20 solved according to the computation model in a similar manner. 
The heat balance for the gas is given by: 

F~c,Cr tJM -T gJi )-h k A k (T gtk -T sk ) = 0 

25 where T 9ik and T Sik are the gas temperature and catalyst 

temperature respectively in channel segment k, c p is the thermal 
capacity of the gas, h k is the heat transfer coefficient in channel 
segment k, and A k is the wall area in channel segment k. The 
heat balance for the catalyst is given by: 

30 

dT 

m n j 
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where m s k is the mass of catalyst in channel segment k, c PtS is 
the thermal capacity of the catalyst material, and -Atf, is the 
reaction heat for reaction j. 

As a specialist in the field will appreciate, the computation model 
indicated above may be modified in many different ways and it is 
also possible to use a type of computation model other than that 
indicated above in order to obtain required information for 
calculating the accumulation actual value A1 , the accumulation 
setpoint value A2 and, where applicable, the emission actual 
value E1. 

Figs. 5-9 depict flowcharts illustrating the method according to 
various embodiments of the present invention. The method steps 
illustrated in Figs. 5-9 are effected continuously when the 
catalyst 4 and the relating injection device are active. As a first 
step, the current state of the catalyst is calculated by means of a 
computation model of the type indicated above. Information 
obtained by the calculations in the computation model is 
thereafter used for calculating the aforesaid accumulation actual 
value A1 and accumulation setpoint value A2. The accumulation 
actual value A1 and accumulation setpoint value A2 are then 
compared and on the basis of this comparison a control signal S2 
for proportioning of reducing agent is generated. As a final step, 
proportioning and injection of reducing agent in the exhaust line 
are effected on the basis of said control signal S2. The amount 
of reducing agent injected is thus controlled on the basis of the 
conformity between the accumulation actual value A1 and the 
accumulation setpoint value A2. 

According to the embodiment illustrated in Fig. 5, an emission 
setpoint value E2 of the type described above is determined, 
whereby this emission setpoint value E2 is taken into account in 
calculating the accumulation setpoint value A2. 
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An emission setpoint value E2 is also determined in the 
embodiment illustrated in Fig 6. In this embodiment, an emission 
actual value E1 of the type indicated above is likewise calculated 
on the basis of information from the computation model. This 
5 emission actual value El is compared with the emission setpoint 
value E2 and this comparison serves as the basis for generating 
a control signal f sp which is then taken into account in calculating 
the accumulation setpoint value A2. The calculation of the 
accumulation setpoint value A2 is thus affected by the conformity 
10 between the emission actual value E1 and the emission setpoint 
value E2. 

In the embodiment illustrated in Fig. 7, a limitation factor f CO nstrain 
is calculated on the basis of information from the computation 

15 model, which limitation factor has a value which depends on an 
estimate of the current risk that the reducing substance content 
of the exhaust gases leaving the catalyst might exceed a 
predetermined limit value. This limitation factor Uonstrain is then 
taken into account in calculating the accumulation setpoint value 

20 A2 in such a way that the accumulation setpoint value A2 
decreases in response to increasing risk that the reducing 
substance content of the exhaust gases leaving the catalyst 
might exceed the predetermined limit value. The limitation factor 
fconstratn is used with advantage as a multiplication factor in 

25 calculating the accumulation setpoint value A2, and to this end it 
is given a value which ranges between 0 and 1 depending on the 
current risk that the reducing substance content of the exhaust 
gases leaving the catalyst might exceed the predetermined limit 
value, whereby the value of the limitation factor is close to 1 

30 when there is no such risk and close to 0 when such risk is 
imminent. In this embodiment, a control signal f sp of the type 
described above is likewise taken into account in calculating the 
accumulation setpoint value A2. 



35 



In the embodiment illustrated in Fig. 8, the accumulation setpoint 
value A2 is calculated by multiplication of two multiplication 
factors, whereby a first multiplication factor takes the form of a 
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calculated accumulation maximum value A max representative of 
the maximum permissible reducing substance accumulation in 
the catalyst under prevailing operating conditions, and a second 
multiplication factor f sp depends on the conformity between the 
emission actual value E1 and the emission setpoint value E2. 
Said second multiplication factor corresponds to the control 
signal f sp indicated above. The accumulation maximum value 
A max is calculated on the basis of information from the 
computation model. Here again it is advantageous to calculate a 
limitation factor f CO nstrain of the type indicated above, whereby this 
limitation factor f constrain is taken into account in calculating the 
accumulation maximum value A max in such a way that the 
accumulation maximum value A max decreases in response to 
increasing risk that the reducing substance content of the 
exhaust gases leaving the catalyst might exceed the 
predetermined limit value. 

Maximum permissible reducing substance accumulation in the 
catalyst means here the greatest reducing substance 
accumulation which, according to calculations performed, can be 
permitted in the catalyst under prevailing operating conditions 
without the release of unconsumed reducing substance to the 
environment exceeding the limit value determined. 

In the embodiment illustrated in Fig. 9, a limit value profile e^ aJt 
is calculated for the ammonia accumulation along the catalyst, 
which is then used for determining the catalyst's capacity for NO x 
conversion under prevailing operating conditions. The NO x 
conversion capacity determined is used in calculating the 
accumulation setpoint value A2, thereby making it possible to act 
upon the accumulation setpoint value A2 in such a way as to 
reduce the risk that the reducing substance content of the 
exhaust gases leaving the catalyst might exceed a predetermined 
limit value. It thus becomes possible to eliminate or at least 
minimise the risk that undesired amounts of unconsumed 
reducing substance might accompany the exhaust gases leaving 
the catalyst and thereby reach the environment, without any 
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limitation factor of the type indicated above having to be used. A 
disadvantage associated with using the previously described 
limitation factor is that inertia and great non-linearity of 
regulation occur in cases where the limitation factor changes 
5 from low to high value or from high to low value. The 

embodiment according to Fig. 9 is not subject to any such 
disadvantage and enables linear regulation with very good 
performance. 

10 In a simplified variant, the method according to the invention may 
be so designed that the accumulation actual value A1 and the 
accumulation setpoint value A2 refer respectively to the current 
and required reducing substance accumulation in the 
computation model's catalyst segment situated nearest to the 

15 inlet end of the catalyst, i.e. nearest to the latter's upstream end. 
In this case, the accumulation actual value A1 can be obtained 
directly from the computation model, and the accumulation 
setpoint value A2 may for example be obtained by the calculation 
algorithm set out below, according to which the catalyst is 

20 divided in its longitudinal direction into a multiplicity of segments 
in a manner corresponding to the computation model indicated 
above. The calculation algorithm used for determining the 
accumulation setpoint value A2 and the accumulation actual 
value A1 is hereinafter called the proportioning algorithm. 

25 

The first thing to determine is a limit value profile for the mole 
fraction of reducing substance through the catalyst: 

30 

where y LiR is the limit value for the mole fraction of reducing 
substance leaving the catalyst, and y N ox,k and y NOx .K are the 
respective mole fractions of NO x in segments k and K (K = the 
last segment of the catalyst). As the mole fraction for NO x 
35 decreases along the catalyst, the mole fraction for reducing 
substance in the limit value profile will likewise decrease along 
the catalyst. The limit value profile for the mole fraction of 
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reducing substance is then used for calculating a limit value 
profile for accumulation of reducing substance: 

RM 1 + K(T SJl )c totJ< -y L Rik 

5 

where K(T s k ) is the equilibrium constant for adsorption of 
reducing substance at the temperature T s in segment k, and c to t,k 
is the total gas concentration. The equilibrium constant is 
determined by the ratio between the respective velocity 
10 constants k a (T Stk ) and k d (T s>k ) for adsorption and desorption of 
reducing substance: 

15 The ratio between complete NO x conversion and current NO x 
conversion is determined by: 

f _ yNp xl o 

'max ~~ 



yNO x .O yNO x ,K 

20 This ratio is used for determining the profile for the mole fraction 
of reducing substance in the case of complete NO x conversion: 

y^k — ytJR + fmax {/nO x M " > 'no x J< ) 

25 This profile is then used for determining a corresponding 
reducing substance accumulation profile: 

gmax _ ^(T s jk ) ' C tot,k ' M 

RM 1 + K(T sM )-c iotJ< -y™ 



30 



Desorption velocities are calculated both for current reducing 
substance accumulation and for the limit value of reducing 
substance accumulation: 
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r^=K(T sM )0 L R , 



5 



Thereafter, the ratio between these velocities is calculated: 




10 



The ratio is used for determining a limitation factor: 



constrain 



1 + (c,% 




.turn 



+ c 



where segments a and b and the parameters cf e ,cf e and cT e 
15 can respectively be selected and trimmed in order to achieve 

optimum functioning. The limitation factor has the characteristics 
of being close to unity when the accumulation of reducing 
substance is small relative to the accumulation limit value. When 
the accumulation is close to the limit value, the value of the 
20 limitation factor decreases. The accumulation setpoint value A2 
is finally determined by: 



25 where f S p is the control signal from the outer regulating circuit 
12b. Here it may be seen that the limitation factor affects the 
accumulation setpoint value A2. When the reducing substance 
accumulation is close to the limit value, the accumulation 
setpoint value A2 will thus decrease until a state of equilibrium 

30 obtains. In this case the accumulation actual value A1 is the 
accumulation of reducing substance in the first segment of the 
catalyst: 




Al = 0, 
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In the proportioning algorithm formulae set out above, R denotes 
the current reducing substance. In cases where the reducing 
substance takes the form of ammonia, R therefore stands for 
5 NH 3 . 

With this simple variant of the proportioning algorithm it is 
possible to achieve rapid regulation of the accumulation of the 
reducing substance in the catalyst. As it is the reducing 

10 substance accumulation in the catalyst's first segment that is 
regulated, the response to changes in proportioning of reducing 
agent is rapid. It is also possible to trim the regulation so that 
the amount of reducing substance leaving the catalyst without 
reacting is constantly below the limit value, since a limitation 

15 factor lowers the accumulation setpoint value A2 when the 

accumulation of reducing substance is close to the limit value. 
The principal disadvantage of this simple proportioning algorithm 
variant is that the internal control state, i.e. the reducing 
substance accumulation in the catalyst, is not directly 

20 representative of the NO x conversion in the catalyst. As the NO x 
conversion is temperature-dependent, more accumulation at 
lower temperatures is required for achieving the same NO x 
conversion. This means that it is the outer regulating circuit 12b 
that has to compensate the accumulation setpoint value A2 in 

25 response to temperature changes. Another problem is that 

temperature changes cause the temperature along the catalyst to 
vary and the reducing substance to be consumed at different 
rates in different parts of the catalyst. As the internal control 
state is only related to the first segment of the catalyst, the 

30 variations across the catalyst have to be compensated by the 
outer regulating circuit 12b. 

In a more refined variant, the method according to the invention 
may be so designed: 
35 - that for each of the segments of the computation model an 
accumulation actual value A k and a conversion value R m ax,k are 
calculated, whereby the accumulation actual value A k is 
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representative of the maximum reducing substance accumulation 
in the segment which is permitted under prevailing operating 
conditions, and the conversion value R ma x,k is representative of 
the expected exhaust gas substance conversion in the segment 
5 when the reducing substance accumulation in the segment 
corresponds to the accumulation value A k , 

- that the conversion values R max , k for the various segments are 
summated, and 

- that the resulting sum is converted to a fictitious value for the 
10 maximum permissible reducing substance accumulation in the 

segment situated nearest to the inlet end of the catalyst, whereby 
this fictitious value constitutes the aforesaid accumulation 
maximum value A max . 

15 The accumulation maximum value A max thus obtained then serves 
as the basis for arriving at the accumulation setpoint value A2 in 
the manner indicated above by multiplication by a multiplication 
factor f sp which depends on the conformity between the emission 
actual value E1 and the emission setpoint value E2. 

20 

A limitation factor f C onstrain,k of the type described above is with 
advantage calculated for each of the segments, whereby this 
limitation factor f C onstrain,k is taken into account in calculating the 
conversion values R max ,k in such a way that the conversion 
25 values decrease in response to increasing risk that the reducing 
substance content of the exhaust gases leaving the catalyst 
might exceed the predetermined limit value. 

A value R K for the current conversion of the exhaust gas 
30 substance in the segment is further calculated for each of the 
segments. A value R to t for the total current conversion of the 
exhaust gas substance in the catalyst is then calculated by 
summating the values R k of the various segments, followed by 
converting the value R toi thus obtained for the total current 
35 conversion of the exhaust gas substance in the catalyst to a 
fictitious value for the current reducing substance accumulation 
in the segment situated nearest to the inlet end of the catalyst, 
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whereby this fictitious value constitutes the accumulation actual 
value A1. 

In this refined variant, the internal control state provides a 
5 representative measure of the NO x conversion through the whole 
catalyst. In this case the proportioning algorithm will be related 
to how the computation model for the catalyst is formulated. The 
following example illustrates how the proportioning algorithm 
might be designed if the reducing substance is ammonia and the 
10 computation model for the catalyst comprises the following 
reactions with relating reaction velocities: 



S + NH 3 -> S-NH 3 






(1) 


15 S-NH 3 -> S + NH 3 


r 2 


— ^2^NH 3 


(2) 


4 S-NH 3 + 4NO + 0 2 ->4S + 4N 2 + 6 H z O 


r 3 


= ^3 C NO^NH 3 


(3) 


4 S-NH 3 + 5 0 2 -4 4 S + 6 H 2 0 + 4 NO 


r 4 


- k 4 c O2 0 NH3 


(4) 



The reaction velocities are given per unit mass of catalyst. 

20 Summating the reaction velocity for reaction 3 (NO x conversion) 
multiplied by the catalyst mass in the respective segments, and 
then dividing the result by the expression for the reaction velocity 
in the first segment multiplied by the catalyst mass but with the 
accumulation of ammonia omitted, provides a measure of the 

25 NO x conversion through the whole catalyst, expressed in 

ammonia accumulation for the first segment in the catalyst, which 
constitutes the current value of the internal control state, i.e. the 
accumulation actual value A1 : 

K 

A\ = <3? = — 

k 3 (T s4 ) c Mtl -y NOXtk -w, 

30 

The setpoint value for the internal control state, i.e. the 
accumulation setpoint value A2, is determined in a similar 
manner: 
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K 



2 f construing ' k 3 ( T s,k ) " <W ' JJVor.* * ' W * 

_ f Jt=l . 

SP,NH 3 JSP' . (T v 



^2 = *^=/, * =J 



where 0^ may be determined in the same way as 0]^ in the 
proportioning algorithm described above in the simplified variant 
of the method according to the invention. The limitation factors 
may be given different values for each segment according to: 



f constraitiyk ■* tune f 

1 + f dJk 



10 where f dik is determined in the same manner as in the 

proportioning algorithm described above in the simplified variant 
of the method according to the invention, and where the 
parameter cf ne may be trimmed differently for each segment. 



15 This refined variant has the same advantages as the simplified 
variant. It is possible to achieve rapid regulation of the internal 
control state, since the latter is closely related to the reducing 
substance accumulation in the first segment of the catalyst. The 
proportioning algorithm can also be trimmed so that the amount 

20 of reducing substance leaving the catalyst without reacting is 

constantly below the limit value, since the limitation factors lower 
the accumulation setpoint value A2 when the reducing substance 
accumulation is close to the limit value. In contrast to the 
simplified variant, the internal control state in the refined 

25 proportioning algorithm is related to the NO x conversion. This 
means that the accumulation setpoint value A2 is automatically 
adjusted when the temperature changes. The fact that the 
internal control state is related to the NO x conversion through the 
whole catalyst means that the proportioning algorithm is also 

30 able to compensate for varying consumption of reducing 
substance in different parts of the catalyst, caused by 
temperature variations along the catalyst. The refined 
proportioning algorithm is thus able to compensate for most of 
the variations which may occur along the catalyst, which means 
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that the outer regulating circuit 12b need only make minor 
adjustments to the accumulation setpoint value A2 in order to 
hold the NO x emission at a desired level. 

A possible implementation of the embodiment according to Fig. 9 
is described in more detail below on the basis that the reducing 
substance is ammonia, whereby the computation model for the 
catalyst comprises the following reactions with relating reaction 
velocities: 

S + NH 3 -» S-NH 3 r, =/c,c NH ^ v 

S-NH 3 -» S + NH 3 r 2 =k 2 0 N „ 3 

4 S-NH 3 + 4 NO + 0 2 -> 4 S + 4 N 2 + 6 H 2 0 r 3 =k 3 c NO 0 NH3 

15 Information from the catalyst model is used for estimating the 
NO x conversion capacity of the catalyst. Surface concentration 
factors q for NO x and ammonia are calculated as follows: 



10 



20 



VnO„Jc ~ 

yNH 3 ,k,s 



where yNOx.k.s and yNOx.k.g are the respective NO x mole fractions 
on the internal wall surface and in the gas flow in channel 
25 segment k, and yNH3,k,s and yNH3,k,g are the respective NH 3 mole 
fractions on the internal wall surface and in the gas flow in 
channel segment k. The denominators are monitored if they are 
nil and in such cases the surface concentration factor is set at 1, 
i.e.: 

30 



W,l =0 = 1 



WO 2005/068797 



28 



PCT/SE2004/001984 



An ammonia concentration factor fNH3.k is calculated as follows: 



where yurea,k,g is the mole fraction of urea in the gas flow in 
channel segment k. 



An iterative procedure is used thereafter for determining the limit 
value profile 0^ M for the ammonia accumulation along the 

10 catalyst. In the first iteration, the derivative for the molar limit 
value fraction of NO x relative to the molar limit value fraction of 
ammonia is set at an initial value: 



15 



20 



dy L N 



NO, 



= -0.5 



In the subsequent iterations, this derivative is calculated as 
follows: 

d0 N Hl fc 3i , (y NOxfi - y NOxA )+ y NOxtl ■ 0 NHiA f 



del 



3 



K <T sA ) • c toUl ■ rj mjA (k 2 (T sA ) + iv k 3 (T sA ) • c MA • y L m ^ • r/^) 



dy L N 



NO, 



dy L NH .. 



3 J 



dy L Nl 



NO, 



<*0 N „ 3 dy L NH3 



dy NOt 
dK H3 



dKn 3 



-1< 



<-1 

dyL 3 
dy L NO , de NH 



d0NH 3 dy N 



<0 



dy L N 



NO x 



d0 L NHa 



d#U 



dy L N , 



>0 



NH, 



25 
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10 



15 



20 



25 



dy L m 



NO x 



dy L Ni 



NH 



i-f-y 



In each iteration the following calculations are also done: 



/ 



L, update 





ay™, . 


i 

•( 2 yi.,NH } + yNo x .t -ym^) 


yNo^o y^j + 




i 

\y/to„i y no j ,k ) 



0 

^ L, update 
yNO x ,Q 



jr L.update < q 

0<f L ' update <- 

^* L, update ^ 



ym x ,o 



yNo xt o yNO x ,K 

y^ro ri o 



yNO xt 0 yNO xt K 



yNo x ,o yNO x ,K 

ym^k -yi t NH 3 + f L * fm^k ~{yNo x ,k ~ym x> K) 

yNO x * - 3>M> a ,0 / Z " 6^,0 "" ^a x ,A ) 



a XH^k 



In the above formulae, y LNH ^ Is the permissible limit value for the 

mole fraction of ammonia after the catalyst, T Stk and c totfk 
respectively are the catalyst temperature and the total gas 
concentration in channel segment k, v 3 is the stoichiometric 
coefficient for ammonia in the NO x reduction reaction, and f L is a 
capacity factor which describes the ratio between NO x conversion 
at the catalyst's capacity limit and current NO x conversion. The 
calculation converges quickly and four iterations are sufficient. 
An iteration factor f ' is used for stabilising the calculation, and 
suitable values for this iteration factor in iterations two to four are 
f 2 =0.2, f 3 =0.1 and f 4 =0. 
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The resulting limit value profile tf^for ammonia accumulation 
along the catalyst can thereafter be used in regulating in order to 
determine the NO x conversion capacity of the catalyst under 
current operating conditions. This capacity is taken into account 
in calculating the accumulation setpoint value A2. Examples of 
conceivable formulae for determining the accumulation actual 
value A1 and the accumulation setpoint value A2 when 
implementing the embodiment according to Fig. 9 appear below. 



In this case the accumulation actual value A1 is calculated from 
the NO x conversion in the first wall segment in all the channel 
segments along the catalyst, converted to an ammonia 
accumulation in the first wall segment of the first channel 
15 segment, according to the following formula: 

K 



* 3 (X.! ) • c tot>1 ■ y NOx>l ■ Tj NO „ k - w, 



where w k is the mass of active catalyst material in channel 
20 segment k. In the formula set out above, the numerator provides 
a measure of the catalyst's current NO x . 

The accumulation setpoint value A2 is obtained by an analogous 
calculation using the limit value profile d^ iJc for the ammonia 

25 accumulation: 



M, JSP MW-WyMW-VMV-". 



30 



In the above formula, the numerator provides a measure of the 
catalyst's NOx conversion capacity under prevailing operating 
conditions. 
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The invention is of course in no way limited to the preferred 
embodiments described above, since a multiplicity of possibilities 
for modifications thereof are likely to be obvious to a specialist in 
the field without having thereby to deviate from the basic concept 
5 of the invention such as it is defined in the attached claims. The 
exhaust system may for example comprise at least one additional 
catalyst connected in series with the reduction catalyst 4, e.g. an 
oxidation catalyst and/or a hydrolysis catalyst upstream from the 
reduction catalyst and/or a slip catalyst downstream from the 
10 reduction catalyst. 



